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Experiment 865 at the Brookhaven AGS obtained 410 K + — > e + ue + e~ and 2679 K + — > [i + ve + e~ 
events including 10% and 19% background. The branching ratios were measured to be (2.48 ± 
O.U(stat.)±O.U(syst.J) x 10" 8 (m ee > 150 MeV) and (7.06±0.16±0.26) x 10~ 8 (m ee > 145 MeV), 
respectively. Results for the decay form factors are presented. 

PACS numbers: 13.20.Eb, 13.40.Ks 



Chiral Perturbation Theory (ChPT) Q has been a suc- 
cessful approach to describing the decays of pseudoscalar 
mesons. In the ChPT program radiative kaon decays 
can serve both as an important test and a source of in- 
put parameters for the theory. While the decay modes 
K + — > e + vj (K e 2 7 ) and /i + ^7 (A^2 7 ) have allowed some 
study of the form factors involved ||, [|| f|, the decays 
K + — > l + ve + e~ (K e 2ee, Kfj,2ee) allow a more detailed 
investigation into the structure of these decays. We re- 
port here on such an investigation from Experiment 865 
at the Brookhaven National Laboratory AGS with a 100- 
fold increase in the number of events in the former mode 
and 150- fold increase in the latter ||. 

The K e 2ee and K^ee decays are assumed to proceed 
via exchange of a W 7+ -boson (l + v) and photon (e + e~). 
The decay amplitude j^, Q includes inner bremsstrahlung 
(IB) corresponding to the tree diagrams in Fig. |l|a,b, 
and structure dependent (SD) radiation (Fig. |]c) P a ~ 
rameterized by vector Fy, axial Fa and R form factors. 
K e 2~/ and K^-y experiments were actually sensitive only 
to \Fy + Fa\- R, which contributes only to decays with 
an e + e~-pair, has not yet been measured. 

Inner bremsstrahlung is unambiguously predicted by 
the K — > Iv amplitude and is proportional to the kaon 
decay constant Fk = 160 MeV. We included the kaon 
electromagnetic form factor in Fig. |l|b (amplitude A4 of 




e + e 



b) c) 

+ ve + e~ decay diagrams 



Ref. 0) in our definition of the IB term. The IB am- 
plitude is negligible in K e 2 ee decay due to electron helic- 
ity suppression, but dominates in K^ee- It contributes 
about 60% of the total K^ee branching ratio for invari- 
ant masses m ee > 145 MeV. An additional 20% comes 
from the IB and SD amplitude interference, which makes 
it possible to measure the signs of all form factors relative 
to F K . 

Generally form factors depend on W 2 and q 2 , where W 
and q are 4-momenta of the Z + z/-pair and of the photon 
(e + e~ pair), respectively. In our analysis we assume the 
dominance of low lying resonances H : 



F V%R - F V,A,R I [(1 - «(1 - W/m 2 )] 



(1) 



where 



= 770 MeV, and in = 



= 892 MeV for 



F v and m = m Kl = 1270 MeV for F A , R. Only the 
constants Fy, Fa, R, which we define in accordance with 
the Particle Data Group (PDG) ||, will be the subject 
of our analysis. The estimated uncertainties in the slope 
of form factors are taken as the model errors. 

We also included in the analysis a hypothetical tensor 
amplitude 



ieG F V u , 
V2 



F T e»q p u v {l + Y 



(2) 



since a possible tensor interaction has been discussed in 
regard to the ir — > ev"f [|| and K — > evir° |hJ experi- 
ments. 

The experimental apparatus was constructed to search 
for the decay K + — > Tt + fi + e~ in flight from an un- 
separated 6 GeV/c K + beam, and has been described 
elsewhere |ll]]. The K e 2 ee and K^2ee data were ob- 
tained in a 1996 run simultaneously with that for a 
measurement of K + — > n + e + e~ (K^ ee ) frEq. The 
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FIG. 2: Missing mass distributions for K^ee and K e 2 ee 
decays. Error bars are data, dashed lines indicate background, 
and solid lines represent simulated distributions. 



trigger for these modes allowed us to preselect events 
with three charged tracks, including an e + e~ pair with 
high invariant mass m ee . Prescaled decays K — > 
Tnr° D {K„ 2 ), fivir° D (K^), evn D {K e3 ), tttt ^ (K„ 3 ) fol- 
lowed by Dalitz decay it a D — > e + e~7 with low m ee were 
used for normalization. 

Off-line the K e 2 ee and K^ee candidates were required 
to have small missing neutrino mass, m u . This missing 
mass was calculated using the measured decay product 
momenta and identities, and the centroid of the incident 
kaon beam momentum determined from K + — > 7r + 7r + 7r~ 
decays. After correction for the decay vertex dependence 
of the momentum, the kaon momentum resolution was 
a p/p — 1.3% and the angular resolution ag x = ag y = 
4 mrad. The m 2 distribution for each decay modes is 
displayed in Fig. |. A cut |to 2 | < 0.016 GeV 2 isolated 
Ki2ee decays. 

Background in both cases was dominated by accidental 
overlap tracks. Events with one of the tracks out of time 
gave model independent samples of the accidental back- 
ground. The samples were normalized by counting events 
with m 2 < —0.03 GeV 2 with normalization uncertainties 
8% for K^ee an d 25% for K e2ee . Other backgrounds and 
processes of interest were simulated using a GEANT3 |1| 
based Monte Carlo. 

The K — > 7ree decay is a potential background both 
for K e 2ee and K^ee decays. To suppress it , each event 
was tested as K ve( .. Events giving effective 7ree mass 
and total momentum close to the beam kaon values were 
removed. 

Cuts on the invariant e + e~ mass > 145 MeV (K^ee) 
and > 150 MeV (K e 2 ee ) removed backgrounds associated 
with large branching ratio processes including a low mass 
e + e _ -pair, e.g. K ^ itTt D . 

After the m ee > 150 MeV cut 410 detected K e 2 ee can- 
didate events remain, including an estimated 35 acciden- 
tal and 5 K^ ee background events. The normalization 
sample was 86k K e3 events and 2.3k events. The 
accuracy of the normalization was estimated to be 4%, 
including the 1.5% error in the K e3 branching ratio, a 



1% trigger efficiency error, 2% for the radiative correc- 
tion, and 3% for reconstruction efficiency. In evaluating 
the normalization factor we assumed that the tt° — > eej 
branching ratio is 1.184%, from the QED calculation Jl4|] 
with a form factor slope 0.032 ± 0.004 |§. 

The scintillation hodoscope embedded in the muon 
stack behind 40 cm of iron Jll| was used for additional 
[ij'K separation in K^ee analysis. The /i/tt detection ef- 
ficiencies in the hodoscope were studied using K„2 and 

events where the ir^ was fully reconstructed, i.e., 
the resulting photon was observed and included in the 
kinematic reconstruction. The average efficiency of de- 
tecting muons with momentum greater than 0.9 GeV/c 
was 90%, while the probability of the misidentifying the 
pion as a muon was about 20%. Because GEANT sim- 
ulation predicts a larger pion misidentification, the cor- 
rection factor 0.85 was applied to the simulated efficiency 
of the pion detection. To reduce the uncertainty of this 
correction, the K^2 decays in the normalization sample 
were additionally suppressed by a cut on < 200 MeV 
(calculated in the kaon center of mass). The final nor- 
malization sample contained 20.5k events including 16.2k 

, 3k if OT 3, 0.8k K^2, and 0.4k accidentals. The nor- 
malization accuracy was estimated to be 5% including 
2.5% from the branching ratio, 1% from the trigger 
efficiency, 2% from radiative corrections, 3% from recon- 
struction efficiency and 2% from efficiency of fx/ir separa- 
tion. The total number of the selected K^ee events with 
m ee > 145 MeV was 2679, including an estimated back- 
ground of 355 accidentals, 126 K^ ee , and 33 K — > mr^n^ 
events. 

In order to fit the data, we have used the logarithmic 
likelihood function 

C = 2 \m l - m + rrij ln(nj/m,)] 

i 



er 



(3) 
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Here, rrii and rii(Fy, Fa, R, Nmc, Na) are the measured 
and calculated numbers of events in i-th bin of the 5- 
dimensional phase space (247 K^ee and 147 K e 2 ee bins). 
We explicitly included the uncertainties, <jmc and a a, 
in the Monte Carlo and accidental normalizations, Nmc 
and Na, respectively in the fit. Nmc and Na were re- 
garded as independent parameters, while N^ c and N^ 
are the expected values from our studies discussed above. 

The consistency between data and simulation was eval- 
uated by significance level (SL) |||, i.e the probability 
that a random repeat of the experiment would observe a 
smaller C, assuming the model is correct. To calculate 
the SL we simulated the expected distribution of C. In 
this calculation we explicitly accounted for possible vari- 
ations of rii due to the finite Monte Carlo and accidental 
statistics. 

Fit results are summarized in Table || In the combined 
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TABLE I: Fit to E865 data. The errors are ±stat. ± syst. ± model. Results for the form factors Fv,Fa, and R (in units of 
1CP 3 ) are given for the fixed Fk = 160 MeV and Ft = 0. Because the measurements of the form factors are correlated their 
linear combinations are also presented. Branching ratios BR (in units of 10 -10 ) were actually measured for m ee > 145 MeV 
(K^ee) and m ee > 150 MeV (-Ke2ee). All other values of BR are extrapolations. 



K 



jj, + ve + e 



K+ 



e ve e 



Combined Fit 



Expected 



Fv 
Fa 
R 



124 ± 
31 ± 
235 ± 



19 ± 13± 4 
21 ± 14± 5 
25 ± 14± 12 



87 ± 30 ± 8 ± 5 
38 ± 29 ± 11 ± 3 
227 ± 20 ± 10 ± 8 



112 ± 

35 ± 
227 ± 



15± 10± 3 
14± 13± 3 
13± 10± 9 



41 ± 6 6 
230 ± 34 



F v + F A 
F v - F A 
R + F v 
R- F v 
R + F A 
R-Fa 



155 ± 
93 ± 
359 ± 
111 ± 
265 ± 
204 ± 



25 ± 21± 

32 ± 17± 

36 ± 20± 

26 ± 18± 

9± 14± 

46 ± 25± 



5 
7 
14 
11 
7 
17 



125 ± 
50 ± 
314 ± 
139 ± 
265 ± 
189 ± 



38 ± 
44 ± 
34 ± 
37 ± 
14 ± 
48 ± 



12 ± 
15 ± 

11 ± 

12 ± 
10 ± 
18 ± 



3 
7 

12 
5 
6 

10 



147 ± 
77 ± 
338 ± 
114 ± 
262 ± 
191 ± 



21± 15± 

20± 19± 

19± 15± 

20± 14± 

6± 9± 

27± 22± 



6 
11 

8 
6 

12 



±|144±9| 
102 ± 74 1 



SL 



11% 



36% 



12% 



F K (MeV) 
F T 



157 ± 
-6 ± 



7 ± 
13 ± 



5± 
8± 



0.3 
1 



(32 2 ± 37 2 ± 25 2 ± 



~2xl/2 



157 ± 

-4 ± 



5± 
7± 



4± 0.2 
7± 0.4 



160 




-^total 

B-Rm cc >0.140 
BRm a e>0.145 
-B-Rm ee >0.150 



(793 ± 
706 ± 
(628 ± 



18± 28± 0.5) 9 
16± 26± 0.4 
14± 23± 0.3) 9 



(1730±g£ ± 90 ± 80) 9 

(291± 16 ± 17 ± 0.7) 9 

(270± 15 ± 16 ± 0.4) 9 

248± 14 ± 14 ± 0.2 



1300 ± 400 / 300± i°o 1 



a Theoretical value (axial anomaly) Fv/Mk = 
6 ChPT to 0(p 4 ) extrapolation from ty — » eu-f 
c R/Mk = {l/S)F K {r K 
d K — > eu-y and K 
e K 



8n 2 F |0. 



, with the experimental value (r K ) = 0.34 ± 0.05 fm Jl5|. 
experimental data Q ^, 0|, corrected with the slopes ([lFof the form factors 
/ivy experimental data 



+ 



* Previous K 
9 Extrapolated value. 



/K 



experimental data |@]. 



fit the likelihood function used was the sum of K^ee and 
K e 2ee likelihood functions (|J) . Expectations for the mea- 
sured values, included in Table |, are based on the previ- 
ous Ki2-y and Knee experiments, ChPT extrapolations of 
the 7T — > evj measurements, and theoretical predictions. 

We considered three main contributions to the system- 
atic error: (i) uncertainty of the detector efficiencies; (ii) 
uncertainty of the normalizations; (in) errors due to the 
finite statistics of the accidental and Monte Carlo sam- 
ples. All three contributions are comparable, and their 
quadratic sums are shown in Table |. 

The model errors in Table [j] correspond to the uncor- 
rected sum of the 30% possible fluctuation of each slope 
of the form factors given in Eq. (|l]) . Uncertainty in the 
dR/dq 2 dominates the model errors of all form factors. 

We include in Table | both measured and extrapo- 
lated branching ratios. The K^2ee total branching ratio 
is expected to be 2.5 x 10~ 5 Q. Because of our cut of 
m ee > 145 MeV, and because K^ee has a large proba- 
bility of events at low m ee due to its being dominated by 
IB, we are unable to determine the total branching ratio 
for that mode. 



To evaluate the sensitivity of our data to the IB term, 
we have made a fit in which Fk in the definition of IB 
amplitude was regarded as a free parameter. The re- 
sulting value of Fk in Table |l| being consistent with the 
expected value, also serves as a check of the normaliza- 
tion based on the decay. The normalization based 
on the K e 3 decay can not be checked in the same way in 
the K e 2ee analysis. However, we can alternatively deter- 
mine that normalization by considering it as an uncon- 
strained parameter in the combined fit. The ratio of this 
normalization to that obtained in K e ^ was found to be 
(0.93 ±0.12)/(1.0± 0.04). 

Using the current algebra relationship || between 
form factor R and kaon charge radius (r K ), R = 
{1/2>)MkFkv" k )i where Mk is the kaon mass, we can 
calculate (r K ) = 0.333 ± 0.027 fm 2 in agreement with 
the direct measurement 0.34 ± 0.05 fm 2 The dif- 

ference between our value and the experimental value 
of the pion charge radius (r 2 ) = 0.439 ± 0.008 fm 2 (f| 
does not agree well with the ChPT prediction to 0{jr): 
(rl) - (r 2 K ) = (1/32tt 2 F 2 ) In (Af|/m 2 ) = 0.036 fm 2 @, 
where F = 92.4 MeV. 
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FIG. 3: Likelihood functions for form factors Fv, Fa, and 
R. 

If we assume that the form factors are constant instead 
of varying according to Eq. (fj]) we obtain for the mean 
values: F v = 0.131, F A = 0.034, R = 0.257. The value 
of C indicates that constant form factors are less likely 
by a factor of 3.6 than those of Eq. ([!]). 

While measurements of the K e 2 ee decay can only de- 
termine the signs of the form factors relative to one an- 
other, due to the interference between the IB and SD 
amplitudes, knowledge of K^ee allows us to establish the 
signs relative to Fk- For the combined fit, the discrimi- 
nation against wrong sign combinations is illustrated in 
Fig. | 

The presentation of the tensor form factor for K e 2 ee 
in Table | underlines that it is F T which is measured in 
this mode. Our analysis of the K^ee decay is, however, 
sensitive to the sign of Ft ■ We did not find any evidence 
for the presence of a tensor term Eq. (|J). Our result 
does not rule out the value F T = -0.0056 ± O.OOf 7 || 
invoked as a possible interpretation [fl9[ of the 7r — > evj 
data from Ref. ||. Assuming that the K e z tensor form 
factor is related to Ft as /t//+ = 3.8Fr |2(|, we find 
/t//+ = — 0.02±0. 04 which strongly disagrees with value 
0.53±g;?g ± O.fO of Ref. g§. A value compatible with 
was also obtained in recent Ki 3 experiments pl| . 

To summarize, we have measured the K e 2 ee and K^ee 
branching ratios: (2A8±0.U(stat.)±0.U(syst.)) x 10~ 8 
(m ee > 150 MeV) and (7.06±0.16±0.26) x 10~ 8 (m ee > 
145 MeV), respectively. For the first time all Ki2 ee form 
factors were unambiguously measured: Fv = 0.112 ± 
0.015±0.010±0.003(WfeL>, F A = 0.035±0.014±0.013± 
0.003, R = 0.227±0.013±0.010±0.009. Our analysis was 
especially sensitive to the sum of axial form factors R + 
F A = 0.262 ± 0.006 ± 0.009 ± 0.006. We did not find any 
inconsistency between K e 2 ee and K^ee form factors. The 
measured inner bremsstrahlung contribution Fk = 157± 
5 ± 4 MeV agrees well with the theoretical expectation. 
No evidence of the tensor amplitude was found: Ft — 
-0.004 ± 0.007 ± 0.007. Our study of the form factors is 
more detailed but consistent with previous K(ir) — > lv^/ 
and K + — > l + ve + e~ experiments. 

We thank J. Bijnens for his FORTRAN code describ- 
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